Photo-transmutation based on resonance γ-ray source
Fu Guang-Yong1, 2, Dang Yong-Le1, 2, Liu Fu-Long1, 2, Wu Di2, He Chuang-Ye2, Wang Nai-Yan1, 2, †
College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China
China Institute of Atomic Energy, Beijing 102413, China

 

† Corresponding author. E-mail: wangny@bnu.edu.cn

Abstract

High intensity γ-ray source can be obtained through resonance reaction induced by protons. In this work, the possibility of using such high intensity MeV-range γ-ray source to transmute nuclear waste is investigated through Mont Carlo simulation. 197Au(γ, n)196Au experiment is performed to obtain the transmutation rate and compared with the simulation result. If the current of the proton beam is 10 mA at the resonance energy of 441 keV, with the γ photons emitted from 7Li(p, γ)8Be, then the corresponding transmutation yield for 129I in 2π direction can reach 9.4×109 per hour. The result is compared with that of LCS γ-ray source.

1. Introduction

Transmutation by neutrons, produced by high flux accelerated protons bombarding a heavy metal target, is the primary method to deal with long-lived fission products (LLFPs) with large neutron capture cross-sections. However, for nuclides, whose neutron capture cross sections are small, this method is not very efficient. In addition, (n, γ) reactions may transmute stable nuclides into long-lived radiotoxic nuclides. For example, 133Cs is stable, of which the (n, γ) transmutation product will be 134Cs followed by 135Cs. They are unstable and radiotoxic with long half-lives of 2.06 years and 2.3 million years, respectively.

Photo-transmutation can be a good supplement to neutron transmutation.[1] There are several methods to obtain γ rays: positron annihilation in flight, bremsstrahlung, laser-Compton scattering (LCS), and nuclear resonance reactions. γ-ray beams obtained by positron annihilation in flight is of low intensity, thus it is not appropriate for transmutations. Transmutations of 129I based on a laser-induced bremsstrahlung (LIB) source have been already performed experimentally in Germany,[24] UK,[5] and Japan.[6] Theoretical investigations have been done on LLFPs such as 126Sn,[7,8] 107Pd,[9] 93Zr,[10] 90Sr,[11] 135Cs,[1214] and 137Ca.[15] An LCS γ-ray source can produce γ-rays in a wide energy range, which can be utilized to transmute LLFPs. The Upgraded HIγS, operating since 2009, can produce up to 3×109 photons/s. The energy of photons can be adjusted in 1 MeV–100 MeV. The energy of electron beam is 0.24 GeV–1.2 GeV and the laser wavelength is . The rate of transmutation reaction of 135Cs by the Upgraded HIγS is evaluated in Ref. [16] to reach 1011 per hour. Studies[17,18] show that the LCS-based photo-transmutation is more efficient than LIB. Many experimental[19,20] and theoretical[14,1618,21] studies have been done using LCS γ-ray source to transmute LLFPs.

In the present paper, γ rays emitted from resonance reactions induced by proton is utilized to transmute LLFPs. Nuclear resonance reaction can produce mono-energetic γ rays, which can be used to transmute LLFPs. To do so, the energy of incident protons is supposed to be hundreds of keV to several MeV. Figure 1 shows high power proton beams worldwide (non exhaustive) crossing this region.[2230] The mean beam current of high power proton beams that have already been constructed can reach 100 mA.

Fig. 1 High current proton beams worldwide (non exhaustive) in MeV region.
2. Proton-induced resonance γ-ray source for transmutation

The excited levels of the compound nucleus are discrete and widely spaced when the bombarding proton energies are low. The so-called Breit–Wigner formula, as shown below, gives a good description for the cross section of a single resonance level,

where σ(Ep) is the (p, γ) reaction cross section of the protons whose energy is Ep in laboratory frame, Ecm,p is the corresponding energy of protons in center-of-mass frame, , E0 is the resonance energy. ma, mA are the mass number of the incident proton and target nucleus. is the statistical factor, IC, Ia, IA are angular momentum of the compound nucleus, the proton and target nucleus, respectively. λp is the de Broglie wavelength of relative motion, ΓT is the full width at half maximum (FWHM) of the compound nucleus level. Γγ is the partial width for the gamma ray transition to the ground state and Γp the partial width for proton emission. Thick target γ photon yields of some nuclides are shown in Table 1.[3133]

Table 1.

γ yields of some narrow resonance reactions.

.

The GDR peak energy of most LLFPs is about 15 MeV, therefore it would be better to use the reaction 7Li 8Be to obtain higher transmutation rates. According to the data listed in Table 1, the yield of 14.8 MeV (33%) and 17.6 MeV (67%) γ rays emitted from 7Li when the energy of incident protons is 441.4 keV is 1.90×10−8 per proton in total.[31,32] The yield of 9.17-MeV γ-ray emitted from 13C 14N is 4.7×10−9.[33]

Different compound nuclides have different energy levels ranging from several hundreds of keV to tens of MeV, which cover most Giant Dipole Resonance (GDR) regions. The yields of γ-radiation produced by proton beams in narrow resonance region follow the equation:

where is the number of resonance energy γ rays produced per second in the solid angle , is the intensity distribution against the incident proton energy, which is a function of the depth of the target l, is the differential cross-section when the proton energy is between E and and the emitted γ rays is in the solid angle to , n1 is the number density of the target atom in the converter. l1 is the thickness of the converter, E1 and E2 are the upper and lower bounds of the incident proton energy.

Next, using the γ-radiation produced above, we obtain the transmutation rate without considering secondary process:

where nt is the transmutation rate, is the number of γ photons produced per second in the solid angle corresponding to the transmutation target, is the attenuation coefficient of the target which is relevant to the γ photon energy E0, l2 is the thickness of the target, is the cross-section of the photo-nuclear reactions when the energy of incident γ is E0, n2 is the number of nuclei to be transmuted. However, this is a complex process. γ photons can produce electrons and positrons through pair production process, these electrons and positrons produce photons again by bremsstrahlung, which can also add to transmutation rate. It is better to track all the particles with a Monte Carlo simulation.

From Eqs. (2) and (3), we can increase the transmutation rate nt in the following ways: 1) higher proton beam intensity Ip; 2) reactions with higher cross-section of proton induced γ-rays; 3) larger solid angle of the transmutation target.

3. 197Au 196Au experiment

First, to ensure the accuracy of the calculation, the experiment of 197Au 196Au was carried out in the 2×1.7-MV tandem accelerator of China Institute of Atomic Energy (CIAE).[33] Two gold targets of the same size are irradiated. Figure 2 shows the transmutation process. The γ photons are produced by 13C 14N reaction when the resonant energy of incident protons is 1.747 MeV. 13C is evaporated onto a ϕ10×2 mm gold target up to a mass thickness of . Bombarded by the 9.17-MeV γ photons, 197Au can be transmuted into 196Au through reaction. However, 196Au is unstable, it decays to 196Pt through electron capture process and to 196Hg through β decay.

Fig. 2 197Au transmutation process.

The first sample is radiated for 6 hours and the second for 5.5 hours. Then, the radioactivity of the irradiated gold samples is measured in a low-background lead chamber with anti-Compton method used to reduce the background radioactivity. The detected counts obey the exponential law given by

where P is the generation rate during irradiation and λ the decay constant of of 196Au. is the irradiation time, is the cooling time, and is the measuring time. By measuring the decay counts Nc, the generation rate P can be derived from Eq. (4), thus the transmutation counts N0 can be obtained from . The spectrum is presented in Fig. 3. Peaks at 355.73 keV (Iγ=87%) and 333.03 keV (Iγ=22.9%) originates from 196Au. The result is shown in Table 2. In Table 2, the second column is the irradiation time, the third column shows the counts of peak at 355.73 keV, and the numbers of (n, γ) in 197Au, calculated from experimental data, with self-absorption of the gold targets taken into consideration, are listed in the fourth column.

Fig. 3 γ spectrum of irradiated 197Au, peaks at 355.73 keV (Iγ=87%) and 333.03 keV (Iγ=22.9%) originates from 196Au.
Table 2.

Transmutation numbers of different samples.

.

The thick target γ photon yield of 13C was measured with an HPGe detector (GXM35P4-70). The maximum yield of 9.17-MeV γ photons is (4.7±0.4)×10−9/proton.[33] As the average current of the proton beam during the 197Au transmutation experiment is , the number of incident protons during 20700-s irradiation time is 1.04×1018, producing (4.9±0.4)×109 γ photons.

In the Geant4 simulation, 4.9×109 γ photons are assumed to be emitted in 4π direction from a point source. The point source is on the center of the front surface of a ϕ10×2 mm gold target, with zero distance. The number of transmutation is (6.6±0.5)×105, which is 37% less than the experimental results but with no deviation in magnitude.

4. Geant4 simulation of photo-transmutation yields

In this section, we simulate photo-transmutation yields of nuclear waste like 129I. The photonuclear cross sections of 129I extracted from Geant4 packages are shown in Fig. 4.

Fig. 4 Photonuclear cross section. The red circle, black square, and blue triangle present the cross section of 129I 128I, 129I 127I, and in total, respectively. The unit 1 b = 10−28 m2.

Photoelectric effect, Compton effect, pair production, and photo-nuclear reactions are the principal mechanisms by which γ rays interact with matter. Figure 5 shows cross section of each mechanism when γ rays interact with 129I.

Fig. 5 129I cross section of different mechanisms, the black square, red circle, blue triangle, and pink inverted triangle show a cross section of the photo-nuclear reaction, photoelectric effect (pho), pair production (conv), and Compton effect (Compt), respectively.

Assuming the intensity of the proton accelerator is 10 mA at 441 keV, according to the yields of 441-keV narrow resonance of reaction 7Li 8Be in Table 1, γ-ray yield of 14.8 MeV and 17.6 MeV can reach 1.2×109 per second, and the γ-ray source, produced by the reaction 7Li 8Be, is assumed to be point shape and isotropic in 4π direction, of which 33% are 14.8 MeV and 67% are 17.6 MeV. Next, the 129I target, with width and height are 2 cm, is placed 5 cm away from the point source. The 129I 128I yields against target depth are shown in Fig. 6, the maximum yield can reach 5.7×10−5 per γ. Taking the γ-ray yield into consideration, there are 6.8×104 129I 128I reactions per second.

Fig. 6 Yields of 129I 128I, the red circle, black square, and blue triangle show transmutations rates of 14.8-MeV γ photons, 17.6-MeV γ photons, and weighted average respectively.

Regarding transmutation of other LLFPs based on 7Li 8Be 441-keV resonance γ-radiation in the same condition as mentioned earlier but the point source with zero distance is on the front surface center of the targets and the targets are as thick as 10 cm. The rates of transmutation reactions are calculated and shown in Table 3.

Table 3.

Transmutation rates of different LLFPs.

.

The transmutation rate of 135Cs is 109 per hour is shown in this chart in total. According to Ref. [16], among LCS facilities constructed around the world, the best transmutation ability is of the Upgraded HIγS, which can reach 1011 per hour, the HIγS is 109 per hour, and the NewSUBARU is 108 per hour.

5. Conclusion

Considering a 10-mA proton accelerator, nuclear waste transmutation driven by 7Li 8Be 441-keV narrow resonance γ-rays have been investigated through Geant4 simulation. The maximum transmutation yield of 129I is 9.4×109 per hour. The transmutation yield of 135Cs is obtained to be the same order as HIγS of 109, which is less than that of the upgraded HIγS by one order of magnitude but greater than that of the NewSUBARU by two orders of magnitude. In the future, we will improve the transmutation rates by improving the current of the proton accelerator, searching for higher yields of γ-radiation in GDR region, designing better target structure, and so on. Moreover, when compared with the good direction of LCS γ source, the γ-radiation emitted from narrow resonance is in directions, which is appropriate for large areas of LLFP radiation.

Acknowledgment

The authors would like to thank Prof. Bing Guo, Prof. Hong-Wei Wang, Prof. Wen Luo, and Prof. Shi-Lun Guo for valuable comments. We appreciate the help from Dr. Yang-Ping Shen on the simulation.

Reference
[1] IAEA technical reports series, 2004, No. 435, Implications of Partitioning and Transmutation in Radioactive Waste Management (Vienna, IAEA) p. 101 https://www.iaea.org/publications/7112/implications-of-partitioning-and-transmutation-in-radioactive-waste-management.
[2] Magill J Schwoerer H Ewald F Galy J Schenkel R Sauerbrey R 2003 Appl. Phys. 77 387
[3] Ewald F Schwoerer H Dusterer S Sauerbrey R Magill J Galy J Schenkel R Karsch S Habs D Witte K 2003 Plasma Phys. Control. Fusion 45 A83
[4] Liesfeld B Amthor K U Ewald F Schwoerer H Magill J Galy J Lander G Sauerbrey R 2004 Appl. Phys. 79 1047
[5] Ledingham K W D Magill J McKenna P Yang J Galy J Schenkel R Rebizant J McCanny T Shimizu S Robson L Singhal R P Wei M S Mangles S P D Nilson P Krushelnick K Clarke R J Norreys P A 2003 J. Phys. D: Appl. Phys. 36 L79
[6] Rahman A K Md L Kuwabara S Kato K Arima H Shigyo N Ishibashi K Nakajima J H K Goto T Uematsu M 2008 Nucl. Sci. Eng. 160 363
[7] Irani E Sadighi S K Zare S Sadighi-Bonabi R 2012 Energy Convers. Manag. 64 466
[8] Wang X L Xu Z Y Luo W Lu H Y Zhu Z C Yan X Q 2017 Phys. Plasmas 24 093105
[9] Irani E Omidvar H Sadighi-Bonabi R 2014 Energy Convers. Manag. 77 558
[10] Sadighi-Bonabi R Irani E Safaie B Imani Kh Silatani M Zare S 2010 Energy Convers. Manag. 51 636
[11] Sadighi S K Sadighi-Bonabi R 2010 Laser Part. Beams 28 269
[12] Takashima R Hasegawa S Nemoto K Kato K 2005 Appl. Phys. Lett. 86 011501
[13] Wang X L Tan Z Y Luo W Zhu Z C Wang X D Song Y M 2016 Laser Part. Beams 34 433
[14] Dazhi L I Imasaki K Aoki M 2002 J. Nucl. Sci. Technol. 39 1247
[15] Sadighi-Bonabi R Kokabee O 2006 Chin. Phys. Lett. 23 1434
[16] Zhu Z C Luo W Li Z C Song Y M Wang X D Wang X L Fan G T 2016 Ann. Nucl. Energy 89 109
[17] Chen J Xu W Wang H Guo W Ma Y Cai X Lu G Xu Y Pan Q Fan G 2009 Nucl. Instrum. Methods Phys. Res. Sect. 599 118
[18] Chen J G Xu W Wang H W Guo W Ma Y G Cai X Z Lu G C Xu Y Pan Q Y Yuan R Y Xu J Q Yan Z Fan G T Shen W Q 2008 Chin. Phys. 32 677
[19] Li D Imasaki K Miyamoto S Amano S Mochizuki T 2005 J. Nucl. Sci. Technol. 42 259
[20] Li D Imasaki K Horikawa K Miyamoto S Amano S Mochizuki T 2009 J. Nucl. Sci. Technol. 46 831
[21] Rehman H Lee J Kim Y 2017 Ann. Nucl. Energy 105 150
[22] Biarrotte, J.L., 2013. Proceedings of SRF, Paris, France http://accelconf.web.cern.ch/AccelConf/SRF2013/papers/moiob01.pdf
[23] https://aries.web.cern.ch/content/iphi
[24] http://soreq.gov.il/mmg/eng/Pages/SARAF-Facility.aspx.
[25] Vernon Smith H Schneider J D Sheffield R 2001 Proceedings of the 2001 Particle Accelerator Conference, Chicago, USA 3296 http://dx.doi.org/10.1109/PAC.2001.988089
[26] Kasatov D Kuznetsov A Makarov A Shchudlo I Sorokin I Taskaev S 2014 JINST 9 12016
[27] Sugimoto M Akagi T Ebisawa T 2018 29th Linear Accelerator Conference, Beijing, China 308 https://doi.org/10.18429/JACoW-LINAC2018-TU2A04
[28] http://www.ifmif.org/
[29] Cho Y S Song W S Choi B H 2001 Proceedings of the 2001 Particle Accelerator Conference, Chicago, USA aa aa 2503 http://dx.doi.org/10.1109/PAC.2001.988298
[30] http://www.d-pace.com/
[31] Fowler W A Laurisen C C 1949 Phys. Rev. 76 314
[32] Hanna S S Meyer-Schutzmeister L 1959 Phys. Rev. 115 986
[33] Dang Y L Liu F L Fu G Y Wu D Guo B He C Y Wang N Y 2019 Chin. Phys. B 28 060706